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It is shown that in a number of scalar potentials with an unstable local maximum at the origin chaotic inflation is followed 
by new inflation if model parameters are appropriately chosen. In this model density fluctuation can have a large-amplitude 
peak on the comoving Hubble scale at the onset of the slow-roll new inflation and can result in formation of appreciable 
amount of primordial black holes on astrophysically interesting mass scales. 



I. INTRODUCTION 



If overdensity of order of unity exists in the hot early universe, a black hole can be formed when the perturbed region 
enters the Hubble radius. While the properties of the primordial black holes (hereafter PBHs) thus produced were a subject 
of extensive study decades ago [0,^ , there were no observational evidence of their existence and only observational constraints 
were obtained against their mass spectrum [^,^. Recently, however, possibilities of their existence have been raised from a 
number of astrophysical and cosmological considerations and there is an increasing interest in it. 

For example, they may be the origin of the massive compact halo objects (MACHOs) which are dark compact objects 
with typical mass ~ IMq and make up about 0(10"^) of the critical density [^]. While the primary MACHO candidates 
are substellar baryonic objects such as brown dwarfs, it is difficult to reconcile such a large amount of these objects with the 
observed mass function of low mass stars |^ and with the infrared observation of dwarf component Q , unless the mass function 
is extrapolated to the lower masses in an extremely peculiar manner of Population HI stars are produced abundantly at the 
relevant mass scale Hence we should consider PBHs seriously as the second option. Furthermore this possibility can be 
experimentally tested by observing gravitational radiation from coalescing black hole MACHO binaries by laser interferometers 

®- 

Another interesting possibility is PBHs with mass M ~ lO'^g which are just evaporating now [llOl . It has been argued that 
high energy phenomenon associated with evaporation can explain origin of a class of gamma-ray burst. For this purpose their 
abundance should be around = 10~* in unit of the critical density . 

One may also consider formation of much heavier black holes with mass M ~ 10* Mq. Such black holes are expected to 
exist in the center of AGNs and quasars and act as an engine of their activity. Although it is generally believed that 
these black holes are formed after recombination with a specially arranged spectrum of density fluctuations jl^, it might 
be interesting to consider the possibility that they were of primordial origin, which might lead to a new scenario of galaxy 
formation. 

In any case, in order to produce PBHs on some specific scale, we must prepare density perturbation whose amplitude has 
a high peak of 0(10"^) on the corresponding scale. It is difficult, however, to realize such a spectral shape in inflationary 
cosmology jl^-^, which is not only indispensable to solve the horizon and the flatness problems but also the most sensible 
way to generate density perturbations [hsl, because usual models predict a scale-invariant spectrum. Given the smallness 
of the observed anisotropy of the cosmic microwave background radiation (CMB) [^, the amplitude of primordial density 
fluctuation turns out to be 0(10"^) on any observable scales in most common models. 

It is possible, nevertheless, to realize a non scale-invariant spectrum by choosing somewhat contrived forms of the inflaton's 
potential. While examples of the various perturbation spectrum realized in different potentials is found in |^ , we must admit 
that few of the models giving non scale-invariant spectrum with a single scalar field has a motivation in sensible particle 
physics. In particular, the model proposed by Ivanov et al. in order to produce significant amount of PBHs employs a 
scalar potential with two breaks and a plateau in between. Several other toy potentials have also been studied in 

One can construct more natural inflation models for PBH formation if one allows multiple scalar fields. Some of them make 
use of primordially isocurvature fiuctuations jisj, others include multiple stages of infiation with each regime governed by 
different field p^-pTf. We must assign appropriate form of coupling of the scalar fields in each model, which may not always 
be easy. 



In the present paper, we propose a new scenario of multiple inflation which, unlike previous double inflation models p4[-p8||, 
contains only one source of inflation. In this model we employ in the Einstein gravity an inflaton scalar field, (j>, with a simple 
potential, V[(j>\, which has an unstable local maximum at the origin, such as a quartic double- well potential. This is the same 
setup as the new inflation scenari o llSll , but chaotic inflation is also possible if the 4> has a sufficiently large amplitude initially. 
In fact Brandenberger and Kung |29| studied the initial distribution of a scalar field with such a potential and concluded that 
chaotic infiation is much more likely than new inflation. Thus we also start with chaotic inflation, but show that new inflation 
is also possible when 4> evolves towards the origin after chaotic inflation if the parameters of the potential is appropriately 
chosen so that the scalar field has the right amount of kinetic energy after chaotic inflation in order to climb up the potential 
hill near to the origin and start slow rollover there. Hence in this model the initial condition for new inflation is realized 
not due to the high-temperature symmetry restoration nor for a topological reason but by dynamical evolution of the 
fleld which has already become sufficiently homogeneous because of the first stage of chaotic infiation. We shall refer this 
succession of infiation simply to chaotic new inflation. 

With an appropriate shape of the potential, density fluctuations generated during new inflation can have larger amplitude 
than those during chaotic infiation. Furthermore, since the power spectrum of fluctuation generated during new inflation 
can be tilted, it can have a peak on the comoving Hubble scale when the inflaton enters the slow-rollover phase during new 
inflation. If the peak amplitude is sufficiently large, it results in formation of PBHs on the horizon mass scale when the 
corresponding comoving scale reenters the Hubble radius during radiation domination. We shall show below that such a 
scenario is indeed possible with a simple shape of the infiaton's potential. 

The rest of the paper is organized as follows. In §11 we summarize basic features of chaotic inflation and new inflation 
scenarios and in §111 we calculate formation probability of PBHs in our model. Then we report the results of numerical 
analysis with various potentials in §IV and §V. Finally §VI is devoted to conclusion. 



II. CHAOTIC INFLATION AND NEW INFLATION 



First let us consider two potentials which have a global minimum at = u < Mpi: 
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The former is the Coleman- Weinberg potential which was employed in the original version of the new inflation scenario 
[ p^ , and the latter is a simple double- well potential. For (f> ^ Mpi > v, the both potentials can be approximated at least 
locally by a simple quartic potential V[(j)] = jcf)'^ with A ~ Aln \4>/v\ for the former case, and the evolution of the universe in 
this regime is practically the same as that in chaotic inflation with the quartic potential [Q. That is, from the slow-rollover 
equations of motion. 
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which are valid until the slow-rollover conditions 
the suffix i implies the initial value at t = t 
N{(j> —> (j>e), is therefore given by 
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the latter approximate equality being valid for (f>^ ^ (f>^. 

The amplitude of linear curvature perturbation, $, generated on the comoving scale, I — 27r/fc, is given by 
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Here tk is the epoch when the wavenumber k satisfied k = a{tk)H{tk) during inflation and Sep = H{tk)/{2-K) is the 
root-mean-square amplitude of fluctuation accumulated during one expansion time around t — tk, and f — 3/5 (2/3) in the 
matter- (radiation-) dominated stage. The above expression is valid until the comoving scale I crosses the Hubble radius 
again. In the case with quartic inflaton potential, it reads 



with Nk^ = N{(f){tk) —> (l>e)- The large-angular-scale anisotropy of CMB due to the Sachs- Wolfe effect is then given by 
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In the usual case with single stage of inflation one typically takes Nke — 60 and determine the value of A using the COBE 
data ST/T ~ 1 x lO"'^ to yield A ~ 2 x 10"^^ In the case with second stage of inflation the value of A^^e corresponding 
to the scale observed by COBE can be significantly smaller and the normalized value of A can be different. Taking A^fee = 20 
with another 40 e-folds during new inflation, for example, we find A ~ 6 x 10~^^. 

After chaotic inflation, the scalar field starts oscillation around = u if w is too large, or it overshoots the symmetric state 
(j) — and approaches the other minimum (j> = —v if v is too small. If v is appropriately chosen at an intermediate value, on 
the other hand, it can spend a long enough time near the origin and then slow-rollover new inflation can set in, which ends 
up with either (j> ~ v or (j) = —v depending on the sign of (f) when classical slow-rollover regime begins. We must resort to a 
precise numerical integration of equations of motion to find out what values of v to take. 

Before doing this, let us review the basic properties of new infiation, assuming that the slow-rollover starts from <^ = <^s at 
the epoch t = ts. Since the shape of the potential around the origin is different between the Coleman- Weinberg potential and 
the double-well potential, we must treat them separately. 

First for the former, the potential near the origin can be approximated by 

^/cwM^Vo-^0^ (11) 

with A — \\n\<j)s/v\ and Vo reads Vo — A«*/16. The classical slow-rollover solution satisfies 
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the latter approximate equality being valid for A'^ ^ 1. The curvature perturbation reads 

*(0 = -\/^^l/, (14) 

where Nkf is defined by Nuf = N{(l){tk) <j>f)- Since A is no larger than ~ 20A ~ 10"^" the resultant curvature perturbation 
is at most of order of 10~* for Nkf < 60. Hence chaotic new infiation with the Coleman- Weinberg potential is not suitable to 
produce large enough density perturbation on small scale to warrant significant formation of primordial black holes. 
Next we consider new infiation with the double- well potential (^, which can be approximated as 

VdwM^ Vo-imV', (15) 

with Vo = jv'^ and = Xv^, but for later convenience we treat as if Vo and were free parameters. The classical slow 
roll-over solution reads 
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Inflation ends ai (j) = (j)f = cv with c being a constant of order of 0.1, so we find 



Ho = ~ e " . (16) 



iV(0-.0,) = ^ln^ = ^ln^, (17) 
and curvature perturbation in the slow-rollover regime reads 

The above classical slow-roll solution (|l^) (and also ( |l^ for the Coleman- Weinberg model) is valid only when classical 
potential force dominates diffusion due to quantum fluctuation. In other words, the classical motion during one expansion 
time, A(f> = \^\H~^ , should be larger than the amplitude of quantum fluctuation, S(j> = H/{2n), accumulated in the same 
period. For the potential ( p^ ) this condition reads 

If <j) stays in the region 101 < <j)q for more than one expansion time, quantum fluctuation will start to dominate the dynamics 
and the universe enters a self- reproduction regime which can last eternally |Q. Then the memory of chaotic inflation would 
be washed away and our double inflation scenario would not work. 

Note that the amplitude of density fluctuation evaluated with the slow-roll approximation exceeds unity, AsR,(tfc) > 1, 
for < (j>q. This means that if (j) entered the slow-roll regime in \ij>\ < 4>q, we would find the peak amplitude of density 
fiuctuation larger than unity. Since our goal is to produce a spectrum of density fiuctuation whose amplitude at the peak is 
at most 0(10"^), we should arrange model parameters so that crosses over the region |0| < (j)q rapidly enough to ensure 
A(tfc) ^ AsR(ffc) in this period and to avoid self-reproduction, or (j> should reverse its direction of motion before reaching 
(f> = 4>q- In either case, as long as we choose the model parameters so that the peak amplitude of fiuctuation is 0(10"^), we 
can automatically avoid the dangerous self-reproduction in the new infiation regime. 

Equation ( p^ implies the spectral index of power-law density fiuctuation is given by 
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Thus the power spectrum can be significantly tilted. The linear perturbation has a peak amplitude 



n = l-^<l, (20) 
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on the comoving scale h = 2n/ka where ks = a{ts)H{ts). It can be large if (jis turns out to be small. Since the coupling 
constant is determined by COBE observation, we have essentially one parameter, v, which determines not only the ratio 
Hq/tti^ but also 4>s. Hence we must await the result of numerical solution of the equations of motion in order to see if this 
model realizes large enough density perturbation on cosmologically interesting scales. 

Before proceeding to the numerical analysis, we next calculate the abundance, if at all, of the PBHs produced due to the 
potentially large fluctuations generated in the new inflation regime with the potential dla). 



III. ABUNDANCE OF PRIMORDIAL BLACK HOLES 



The initial fraction of the PBHs with mass M, l3{A4), produced during radiation dominated era has been estimated by 
Carr H as 



/?(M)^5(M)exp(-^^^), (22) 
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assuming the density fluctuation is Gaussian distributed with the root-mean-square at the horizon crossing given by 5{M). 
Since PBHs are formed at high density peaks with (5 > 1/3, however, possible non-Gaussian effect on the tail may significantly 
affect the estimation of their abundance, as was first pointed out by Bullock and Primack |Q . Here we estimate the abundance 
of PBHs produced in our model by calculating the probability distribution function (PDF) of the coarse-grained field in terms 
of the stochastic inflation approach [^, following more recent analysis by Ivanov jj^ . 

Since our model predicts a relatively sharp peak at the scale Is corresponding to the onset of slow-rollover we estimate the 
abundance of PBHs on this particular scale with the potential 



The scalar field coarse-grained over a super-Hubble scale, cj>{t,-x.) satisfies the Langevin equation 

3Ho^ + V'M =g(t,x), 
with the temporal correlation function of the stochastic force, g(t,x), given by 

{g(t,x)g(t',x)) = ^5(t-f'). 
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Since the potential is a quadratic function, the PDF of 0, Yl[(j) = <j),t] remains Gaussian as long as the initial distribution 
is also Gaussian. For example, starting with the PDF H[0,ts] = S{(j) — 4>s) and solving the Fokker-Planck equation 



with the boundary condition H[i 
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This does not imply, however, that the resultant abundance of PBHs are given in the form (g^) , because of the nonlinear 
dependence of the metric perturbations. Following Ivanov |p6|, we write the coarse-grained metric in the quasi-isotropic form 



ds^ = —dt^ + {4>{'t; x)) dx^. 



(30) 



where the scale factor a now depends on the coarse-grained spatial coordinate as well, and quantify the metric perturbation 
in terms of 



a{t,yi) 



(31) 



which is frozen in the super-Hubble scales. In the limit h <^ 1, h reduces to the gauge- invariant growing adiabatic perturbation 

We are interested in the statistical distribution of metric perturbation on the comoving scales that leave the Hubble radius 

in the period between t — and r ~ H^^ when the classical solution has rolled down to (j>ci{T — H^^^) — (j>s exp ^^^^ = 0si. 

While the amplitude of the metric perturbation is determined by the duration of inflation in each region, it is not straight 
forward to extract perturbation on these particular scales in the present approach in which nonlinear effects are at least 
partially taken into account. In this situation, in order to obtain the PDF for the duration of inflation in a specific regime, 
Ivanov [M proposed to solve the Fokker-Planck equation (U&) with the initial condition n[0, 0] = S{(f>~(f>s) and the absorbtive 
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boundary condition at (f> = (f>ci{T = Hq and to identify the probability current J in ( |26[ ) with the desired PDF, from which 
he finds the PDF of h, P[h], as 
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in the case with a linear potential, where N is the e-folding number of expansion while the classical solution allows the 
e- folding number of Nd, to which we put unity here, h and TV are related by 



h = exp{N-Nd) 



(33) 



Since our potential is more complicated, it is cumbersome to adopt the same procedure. So we instead use the solution 
( ^t\ ) and obtain the PDF of the duration of inflation while (f> evolves from 0^ to (f>si, r, simply assuming that (f> satisfies the 
classical equation of motion at the moment — (f)si. Then we find the PDF of h as 



P[h = h] = 



mi 1 + h 
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Before proceeding to the calculation of PBH abundance, we mention the relation between this approach and Ivanov's result 
21). If we follow the same procedure as above for a linear potential we obtain explicitly 
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Thus the only difference to ( p^ ) is the absence of a prefactor Nd/N. One may also be tempted to define the PDF for f from 
the probability current J using a solution analogous to ( p^ ) but with a linear potential. In this case one would find 
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Again we only have an extra prefactor of order of unity. Since the deviation between A'^ and Nd is not expected to be too 
large in our case, unlike in the situation of eternally self- reproducing inflationary universe |Q, all the above three approach 
will yield essentially the same result as far as the rare PBH formation is concerned. This also implies that our semi-nonlinear 
approach suffices to the present problem. 

We now proceed to evaluation of the fraction of PBHs produced using (M) , which is explicitly written as 
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Since 4> is slowly rolling at 



one can also write (? in terms of A as 
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In the limit of small h it properly reduces to Gaussian: 
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If we further take the limit 7 



^27rc2(l - e-27) 
0, it reduces to a more familiar form 
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From (H) the probability that h exceeds a threshold value hth is estimated as 

/•oo 

P[h>hth]= / P[h]dh 
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which is insensitive to the upper end of integral since the dominant contribution comes from h ~ hth- The criterion for black 
hole formation has been numerically investigated by Nadegin, Novikov, and Polnarev and by Biknell and Henriksen [ p9| . 
Although it depends on the shape of the perturbed region, the generic value of the threshold reads hth > 0.75 — 0.9. We take 
the black hole threshold as = 0.75. 

Putting hth = hbh we can identify (^) with the fraction of primordial black holes, 

I3{M,) ^ P[h > hih], (43) 

where the typical black hole mass, Ms, is equal to the horizon mass when the comoving scale h reenters the Hubble radius. 
Note that since our model predicts density fluctuation which is highly peaked on the comoving scale h the resultant mass 
function of PBHs are sharply peaked at the mass around Ms- 



IV. NUMERICAL ANALYSIS OF CHAOTIC NEW INFLATION WITH THE SIMPLEST POTENTIALS 

We now proceed to the numerical analysis of the equations of motion for the homogeneous part; 

4>{t) + 3H<j, + V'[(l}]=0, (44) 

with either V[(j}] = Vcw[<?!'] or V[(t>] — Vuw[4>]- One of the nice features of chaotic inflation is that no matter what initial 
configuration is adopted, the system rapidly approaches the slow-rollover solution with exponential accuracy once inflation 
sets in with sufficiently large initial value of (j). We stress in this sense that the dynamics of the second inflationary phase 
is independent of the choice of initial condition but is only sensitive to the model parameters of the potential, and unlike in 
the original version of new inflation fine-tuning of initial condition is not necessary. For our purpose, it is sufficient to take 
(j) > 3.5Mpi -UIpi. 

First, for completeness we solve equations (^) and ( psj ) for new inflation with the Coleman- Weinberg potential (|l|). We 
find that if 0.2201Mp; < v < 0.2259Mpi at least ten e-folds inflationary (accelerated) expansion with \H\ < is realized 
independent of the value of A. For 0.2223Afpi < v < 0.2239Afp( we have more than 60 e-folds new inflation and in this case 
the chaotic inflation regime would be inflated beyond our current Hubble radius. In any case, as we mentioned in §n we 
cannot obtain large-amplitude fluctuations in the observable scale in this model. 

Next we report the case of the double- well potential (^. In this case we flnd, independent of the value of A, that the 
fleld settles down to the positive minimum if v > 0.16286751A/pi and it overshoots the origin to a negative value for 
V < 0.16286750A/pi = Vcr- If v is much smaller than Vcr, 4> will travel between positive and negative values several times 
before settling down to one of the minima, but we do not consider this possibility here. 

For V ~ Vcr, we find the vacuum energy at the origin induces inflationary expansion which lasts only about five e-folds. 
This is because the curvature of the potential at the origin is so large that the universe cannot stay in the slow-roll phase as 
is seen from the fact that the ratio 




is substantially larger than unity. As a result the magnitude of density fluctuations produced in this regime remains small: 
A = C(10~^) with COBE-normalized self coupling A = 3 x 10"^^. 
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We thus find neither the Coleman- Weinberg potential (|l|) nor the double-well potential (g) lead to formation of large- 
amplitude density fluctuations on currently observable scales in the chaotic new inflation scenario. But the reasons of their 
failures are opposite; with the Coleman- Weinberg potential, the field moves too slowly and the scales with large-amplitude 
fiuctuations are inflated away, while with the double-well potential, </> moves so rapidly that A remains too small. 



V. PBH FORMATION IN CHAOTIC NEW INFLATION 



The above observation naturally leads us to consider a different class of potential including another free parameter: 



Vm 



2,2, 

= --m (j) + 



(47) 



that is, typical one-loop effective potential with nonvanishing mass term at the origin. This type of the potential with a 
positive mass-squared at the origin was employed in the original inflation scenario but for our purpose we adopt a 

negative mass term. 

The potential ( ]47| ) has four parameters, but one of them, Vo is flxed from the requirement that the vacuum energy density 
vanishes at the potential minimum (j) = '^4>m- While we numerically solve the equation l^'[(?!>m] = and obtain the numerical 
value of 1^0 so that l^[(j!)m] = 0, we can also calculate them perturbatively in the case ^ , as 
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m 



A 4 1 2 2 



(48) 
(49) 



Another parameter, say A, can be fixed from the amplitude of large-scale CMB fiuctuations using the CORE data as before. 
Hence we are essentially left with two free parameters, v and m. While v mainly controls the speed of (j> around the origin and 
its fate. I.e., to which minimum it falls, and m mainly governs the duration of new inflation, the entire dynamics is determined 
by a complicated interplay of the three parameters. For example, we cannot determine A until we calculate the duration of 
new inflation which also depends on A itself for flxed values of m and v. Hence we must numerically solve the equations of 
motion iteratively to find out appropriate values of parameters to produce PBHs at the right scale with the right amount. 
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FIG. 1. Evolution of the inflaton in chaotic new inflation with A = 3 x 10 
Time and (j) are displayed in unit of the Planck time and Mpi , respectively. 



X 10~^ Mpi, and v = 0.2138436Mp; . 
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Let us now consider a specific example of formation of MACHO-PBHs. For this purpose we must realize a peak with 
/3 ~ 10"^° on the comovine Hubble scale at A'' = 35. After some iterative trials we have chosen A = 3 x 10 and 
m — 6 X 10~^ Mpi, and then solved the equation of motion for various values of v. In this choice of A and m we find new 
inflation lasts for more-than ten e-folds expansion if we take v in the range i; — 0.2131Mp; — 0.2147Mp;. Hence we do not 
need much fine-tuning of the model parameters to realize a new inflationary stage itself. We also find that (p settles down to 
(j)m classically if n > 0.213843638Mpi and to —(f)m if v < 0.213843631Mp; . If, on the other hand, we choose v in the range 
0.213843632A'/p; < v < 0.213843637A/p; , (j) spends more than one expansion time in the region < <j)q = 3.4 x 10"**Mpi 
and the universe would enter the self-reproduction regime. 



Fig. 2 




FIG. 2. Evolution of the scale factor with the same parameters. 



Figures 1 and 2 depict evolution of the scale factor and the infiaton (p, respectively, for the case v = 0.21384360Mp; = vm 
with the initial condition Ui — 1 ai tfji ~ 3.5Afp;. The chaotic infiation ends at <^ ~ 0.89Afpi and new infiationary expansion 
sets in at ~ 0.082Afp; but the slow roll-over phase starts only at (fia = —4.03 x 10~^Afp; = cf)M- In this case is found to 
stay in the region |<ji| < <j)q for only about 0.1 expansion time. The linear perturbation A has the right amplitude on large 
scales to meet the COBE observation, and it has a peak on the comoving Hubble scale ai tf) — 4>s- We find N[(j>s <j}f) = 35, 
which is the right scale for MACHO-PBHs. 

In this case we find 7 — 0.300590 and the abundance of the PBHs at formation reads 

/3 = 2.29c exp(-0.0196994c-^ ) , (50) 

with = 0.75. For (j)^ ~ <f)M we find the peak abundance of PBH, fi — (S x 10"^" at the mass scale M ~ IMc^. Using (ti£l) 
we can also write it as 

/3 = 0.888A«exp(-0.131072A7^). (51) 

One can also obtain an approximate shape of the mass spectrum of the PBHs using (^) with Aj, replaced by A at different 
epoch corresponding to different black hole mass. More specifically the mass of black holes, M, and their initial fraction, 
P{M), can be written by an implicit function of as 



^ In obtaining ( pO| ) we have started with Yl[(j>,t = ts] = 5{(f) — (j)s) in order to extract information on a specific mass scale. 
This also corresponds to treating that </> evolves along the classical trajectory until (j) = (j)s, which is correct only as an average. 
In fact, due to quantum fluctuations generated during chaotic inflation and the early stage of new inflation, 0s itself takes 
different values in different domains. We have confirmed, however, that (ps shifts only about ±4_?/o/(27r) even if we consider 
their effects. Since we find \<j}M\ 3> 4_H'o/(27r), it does not induce significant fiuctuations in /3. For the same reason we are free 
from the domain wall problem which could be present if (j) had a large fiuctuation and different regions fell different minima. 
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M = 7<'exp(2 [N{(t> cj>f) - 35]) Mq, 
P{M) ~ 0.888A(<;/>) exp (-0.131072A"^ (<^)) , 



(52) 
(53) 



where K is a factor of order of unity which depends on the expansion law of the post-inflationary universe but we put K = 1 
for simplicity below. 

Figure 3 depicts the mass spectrum of black holes obtained from (^2|) and (|53[). Thus the PBH abundance is sharply 
peaked. Note, however, that the shape of the large-mass tail is not exactly correct which corresponds to the regime where 
slow-roll solution is invalid. Nonetheless this figure correctly describes the location of the peak. 
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FIG. 3. Expected mass spectrum of PBHs in chaotic new inflation with the same model parameters. Mass is displayed in 
unit of the solar mass. 



Apart from the effects of infiaton's detailed dynamics, however, we must say that the above spectrum is only qualitatively 
correct, because we have chosen a specific value of the threshold, hi,h ~ 0.75, and assumed the black hole mass is equal to the 
horizon mass when the perturbed region reentered the Hubble radius (see also po|pl| ]). In order to improve the calculation 
of the mass spectrum we must calculate the probability distribution functional of the configuration of the perturbed region 
and then calculate the final mass of the black hole, if formed, for each configuration, which is beyond the scope of the present 
analysis. Nevertheless from the above analysis we can convince ourselves that we can produce large enough amplitude of 
curvature fluctuation on a desired scale in the chaotic new inflation scenario. 

Within the limit of our predictability of the mass spectrum, we can also apply our model for the formation of PBHs 
with different masses and abundance. For example, we may produce PBHs with M ~ 10* which may act as a central 
engine of AGNs with the current density, say, n ~ 10~^Mpc~^ corresponding to /? ~ 10^^^ at formation. From ( p^ ) we find 
M = 10* M0 corresponds to N{(f)s (jif) = 44 and the desired spectrum is realized for A = 3.7 x 10~^^, m = 6.1 x 10~*A'/p;, 
and V = 0.21532324Mp, under the COBE normalization. 

Another interesting possibility is to produce a tiny amount of PBHs which are evaporating right now, with the initial mass 
M ~ lO^^g. With the current abundance Q, ~ 10~* or /3 ~ 10~^^ at formation, they may explain a class of gamma-ray bursts. 
In this case we should have only a short period of slow-roll new inflation, N{(j>s (f)f) = 14. We find /3 = 2 x 10~^^ at the 
right mass scale if we choose A = 3 X 10"", m = 5 x lO'^Mpi, and v = 0.16557604828Mpi . In this case a relatively large 
value of m is required in order to keep new inflation short, and we cannot necessarily rely on the stochastic inflation method 
which is valid only if |m^| < Hq. This does not mean we cannot generate large enough density fluctuation on the relevant 
scale. The only problem is we do not have a reliable method to calculate magnitude of fluctuation or black hole abundance 
accurately in such situations. 
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VI. CONCLUSION 



In the present paper we proposed a new scenario of double inflation which contains only one inflation-driving scalar field, 
that is, we pointed out that a scalar potential which has an unstable local maximum at the origin can not only realize new and 
chaotic inflation separately but also accommodate both sequentially if its model parameters are appropriately chosen with 
natural initial conditions as employed in the chaotic inflation scenario. We have further shown that the spectrum of density 
fluctuation in this model can have a large-amplitude peak on the comoving Hubble scale at the onset of the slow-roll regime 
of new inflation and that this can be applied to formation of PBHs on a specific mass scale. This feature of the spectrum is 
realized naturally compared with other models with a single scalar field pi| , p2| because the scalar field is not slowly rolling 
at the onset of new inflation. On the other hand, we must specify the values of model parameters with many digits in order 
to produce an appropriate amount of PBHs on a desired scale. This feature, however, is more or less common to all the 
other models attempting to account for formation of PBHs in inflationary cosmology (see e.g. |^), because both the peak 
amplitude of the fluctuations and its location must be specifled with high accuracy due to the exponential dependence of the 
black hole abundance and its mass on the model parameters. 

On the other hand, one could in principle claim that observation of PBHs can serve as a strong tool to determine the 
parameters in the inflaton's dynamics. Unfortunately, however, our ignorance of the detailed condition for PBH formation, 
such as the precise threshold amplitude of fluctuation as a functional of the shape of perturbed region, makes it impossible to 
link the mass spectrum of PBHs with the shape of the inflaton's potential precisely. In the present paper we have calculated 
the values of the model parameters rather precisely under the universal assumption of hth ~ 0.75. In fact, however, the 
values of the parameters would totally change had we chosen a different threshold. Hence the precise numbers we have quoted 
do not have much significance, but the number of digits simply indicates the sensitivity of the mass spectrum to the model 
parameters. 
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